Random-time processes governed by differential 
equations of fractional distributed order 



L.Beghin* 



Abstract 

We analyze here different types of fractional differential equations, under the assumption that 
their fractional order v 6 (0, 1] is random with probability density n(v). We start by considering 
the fractional extension of the recursive equation governing the homogeneous Poisson process 
N(t), t > 0. We prove that, for a particular (discrete) choice of n(v), it leads to a process with 
random time, defined as A^(7~v, v,(r)), t > 0. The distribution of the random time argument 71,, 
can be expressed, for any fixed t, in terms of convolutions of stable-laws. The new process 
N(Tyj V2) is itself a renewal and can be shown to be a Cox process. Moreover we prove that 
the survival probability of NCTy^ „), as well as its probability generating function, are solution to 
the so-called fractional relaxation equation of distributed order (see ri61). 

In view of the previous results it is natural to consider diffusion-type fractional equations of 
distributed order. We present here an approach to their solutions in terms of composition of the 
Brownian motion B(t), t > with the random time 7~v[ . We thus provide an alternative to the 
constructions presented in Mainardi and Pagnini 1191 and in Chechkin et al. (6), at least in the 
double-order case. 

Key words: Fractional differential equations of distributed order; Stable laws; Generalized 
Mittag-LefHer functions; Processes with random time; Renewal process; Cox process. 
AMS classification: 60K05; 33E12; 26A33. 

1 Introduction 

In the last decade an increasing attention has been drawn to fractional extensions of the Poisson 
process: see, among the others, [261, [10|, [13|, [29], lITSI , ifTTl . flSI, |5|. In particular, the analysis 
carried out by Beghin and Orsingher ||2] starts from the generalization of the equation governing the 
Poisson process, where the time-derivative is substituted by the fractional derivative (in the Caputo 
sense) of order v e (0, 1]: 

^-A{pt-pk-x), k>0, (1.1) 



with initial conditions 



pm-il \l\ (1.2) 



and p-i(t) = 0. The solution to this equation has been expressed as the density of the random-time 
process called Fractional Poisson process (FPP) and defined as 

Ny(t)^N(T2ym t>0. (1.3) 

Here denotes the standard homogeneous Poisson process with rate parameter A > 0, while 
7~2v(f)> f > is a random process (independent from A^) with density given by the folded solution 
to the fractional diffusion equation 
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Later, in [3 J, the distribution of the FPP A/V has been expressed as 



pl(t) = (An'ElX^(-An, k>0,t>0, (1.5) 
in terms of the so-called Generalized Mittag-Leffler (GML) function, which is defined as 

°° (y) ■ 

Kb^^) = ..T., "^P^y^ Re(a),Re{J3),Reiy) > 0, (1.6) 

where (y)j - y{y+ \)...{y + j - 1) (for j - 1,2,..., and y 0) is the Pochammer symbol and (y)o - 1 
(see ifTTI . p.45). Moreover a higher-order generalization of the previous results has been obtained in 
O by introducing "higher-order fractional derivatives" in (11.11 ) and analyzing the following equation 

^ AM-nr& + -A.. M-'^--npu-Pk-i\ k>0, (1.7) 



dt"" \lj dt''"-^> \n-ll dt 

where v e (0, 1) , subject to the initial conditions 



PkiO) = -i J , forO<v<l (1.8) 



= = 1,...,«- 1, k>0, fori<v<l 

f=0 « 



and p-\{t) - 0. The solution to (11.7b was given by the following finite sum of GML functions: 



The corresponding process was proved to be a renewal, linked to Nyit), f > 0, by the following 
relationship 

Wit) = p:,{t) + p:,,,(t) + ... + p:,,„_,(t), t > o. 

Thus it can be interpreted as a FPP which "records" only the ^-th order events and disregards the 
other ones (for an application to the theory of random motions at finite velocity, see |4 |). 

We will introduce here the assumption that the fractional order v of the derivative appearing in 
equation dl.lb is itself random, with distribution «(v), v € (0, 1]: i.e. 

-^n(v)dv = -Aipk - pk-il k>0,ve(0,l]. (1.10) 
dt" 

More precisely, we will concentrate on the case of a double-order discrete distribution of v, i.e. 

n(v) — ni6(v — vi) + n26(v — V2), < vi < V2 < 1, (1-1 1) 

for «i,n2 ^ and such that rti +112 = 1. The assumption ( 11.1 lb has been already considered in the 
context of fractional relaxation (see lfT6l ). as well as for fractional kinetic equations and, in the last 
case, it leads to the so-called diffusion with retardation (see 1 19 |). As we will see in the next section 
this assumption on v produces a form of the solution which is much more complicated than (11.5b and 
( 11.9b . since it involves infinite sums of GML functions. Nevertheless the renewal property is still 
valid and a subordinating relationship similar to (11.3b holds for the corresponding process, which 
can be defined as 

N(fy,,y,(t)), t>0. (1.12) 

In this case the random time Ty,^v2 is represented by a process whose transition density can be ex- 
pressed either as an infinite sum of Wright functions or by convolutions of stable laws. 
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In section 3 we will investigate the relationship between the previous results and the diffusion 
equation of fractional distributed order 



-^n{v)dv ^ X e R, f > 0, v(x, 0) = (1.13) 
at" dx^ 

for < V < 1 . Equations like (11.131 1 have been already studied in ||6] and lfT9l in connection with 
the kinetic description of anomalous diffusions. It has been proved by Chechkin et al. [6| that the 
solution v(x, f), X e R, f > 0, is a probability density function and that the corresponding process is 
subordinated to the Brownian motion via the following relationship 

v(x,f)= I — ^=G{u,t)du. (1.14) 
Jo 



In (11.141) the function G is defined by its Laplace transform 



X{G(m,0;?7) = " e-^io '"''"^^''r (1.15) 

In the special case of double-order fractional derivative in ( 11.13b these authors focus on the behavior 
of the second moment of v(x, t), which suggests that the process can be interpreted as a "diffusion 
with retardation", in this case. Moreover, under assumption dl.l lb . equation ( 11.131) can be seen as 
a particular case (for y = 2) of the equation ( 12.35b below, which is analyzed in 127). In this paper 
only the Fourier transform of the solution is given in explicit form, in terms of infinite sums of 
generalized Mittag-Leffler (GML) functions. Finally the solution to ( 11.13b has been analytically 
expressed in terms of generalized Wright functions by fl9]. 

We prove here that the solution to equation (11.13b with the assumption ( 11.1 lb . i.e. 

m-— + n2^— = T-i", xe R,f > 0, v(x,0) = ^(x), ni,n2 > 0, (1.16) 
ofi of^- dx^ 

for < vi < V2 < 1, coincides with the density of the random-time process 

B(fy,,y,(t)), t>0, (1.17) 

where B is the standard Brownian motion and the time argument 7"v, _v, is the same as in (11.12b . 
thus writing in explicit form the density G in (11.14b as an infinite sum of Wright functions or by 
convolutions of stable laws. _ 

Finally we note that the density of the random time 7"vi _vt appearing in the processes (11.12b and 
( 11.17b coincides with the solution to the equation ( 11.13b . when a difl'erent hypothesis on the density 
n(v) is assumed, i.e. 

n(v) - nl6(v-2vi) + nl6(v -2v2) + 2nin26(v - (vi + V2)), < vi < V2 < 1, (1-18) 

for «i, n2 ^ and such that n\ + n2 - I ■ Therefore the equation governing the process Tv^,v^{t), f > 
turns out to be 



"''dF^"^di^j ^dx^' -^eR,f>0, ni,«2>0, (1.19) 

for < vi < V2 < 1, with the usual initial conditions and v,(x, 0) = 0, in addition. 

Equations (11.16b and ( 11. 19b are proved to govern deeply different processes: while the former 
is linked, for any value of vi,V2 to a diffusion with retardation (see also [6] and fT]), the same is 
not true for the second equation, which, depending on the value of the random indexes, produces a 
subdiffusion or a superdiffusion. 
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2 The recursive equation of distributed order 
2.1 The double-order fractional case 

We begin by considering the following fractional recursive differential equation 



Jo 



Jo df 
where, by assumption, 



-n(y)dv^-A{pk-pk-x), k>0, (2.1) 



n(v)>0, r n(v)t/v=l, v € (0, 1] (2.2) 
Jo 

subject to the initial conditions 

pm = { ; ll'l , (2.3) 

with p-\(f) - 0. We apply in (12.11 1 the definition of fractional derivative in the sense of Caputo, that 
is, for m e N, 

d^ , ^ { r,^ . C 4^u(s)ds, ioxm-\<v<m 

u(t) - { r(m-y) Jo (r-J)i+- ds" v / , ^2.4) 

dt" \ jp^uO), for V = m 

(see, for example, ifTTl . p.92). As a special case, for n(v) = (5(v - v), and a particular value of 
V G (0, 1) , equation (12.11 1 reduces to (II. Il l, which governs the so-called FPP Nv{t), f > (see, for 
details, fT5\, fTT\ and fTSl). 

In order to get an analytic expression for the solution to (12. Il l, we adopt here the following 
particular form for the density of the fractional order v: 

n(v) — ni6(v — vi) + n26(v — V2), < vi < V2 < 1, (2.5) 

for ni,n2 ^ and such that n\ + n2 - I (conditions (12.2b are trivially fulfilled). The density (12.51 1 
has been already used by lfT9l and 16), in the analysis of the so-called double-order time-fractional 
diffusion equation, and corresponds to the case of a subdifFusion with retardation (see next section 
for details). Moreover, it was applied in [16] in the context of fractional relaxation with distributed 
order 

Under assumption (12.5b . equation (12.1b becomes 

d^'pk d^^pk 

n\—r—+n2—r—^-A{pk-pk-\), k>0. (2.6) 



By taking the Laplace transform of (12.6b we get the following first result 



Theorem 2.1 The Laplace transform of the solution to equation ( 12.61 ), under conditions ( 12.31 ), is 
given by 

for any k > 0. 

Proof Formula (12.7b can be easily obtained by applying to ( 12.6b the expression for the Laplace 
transform of the Caputo derivative, i.e. 



m-l jf 

//"X MO; 7?) -^77""'-' — M(f) 



/=0 
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where m - [vj + 1, which yields, for ^ > 1, 



n,7i^'£ [K; ri} + mif^-L {/^; r^] = -A [C {/^; r^]-L[p^,_,■ r^]\ . (2.9) 
By recursively using (12. 9t we get 

^(^^-^H U».,',„,,. f^{/^o;.}.^^i- (2.10) 

For ^ = 0, we get, instead: 

-i- Pn 77 = . (2.11) 

which, together with (12.10b . gives ( I2.7l i. □ 

The Laplace transform (12.7) can be compared with formula (4.8) of 11211 . where a Poisson process 
time-chenged by an arbitrary subordinator is considered. 

We can not use a direct method in order to invert analytically the Laplace transform (I2.7l i. Indeed 
an explicit inversion formula is available only for k - 0, while for ^ > the presence of the power 
k + I makes the analytic inversion too complicated. For A; = 0, we can apply the well-known 
expression of the Laplace transform of the GML function defined in (11.6b (see ifTTI . p. 47), i.e. 

X{f-'4,(-/^);,} = ^^^, (2.12) 

(where Re(y8) > 0, Rs{y) > 0, Re(^) > and i > |w|^) and the resulting formulae (26) and (27) of 
EH. 

Therefore we get 

K(t) (2-13) 

under condition \n\Tf^ I{n2jf- + X)\< I (which is fulfilled, for V2 > v\, A> 0). 

For A; > 0, we adopt an approach similar to those used in lH], ||2], ll22l . ||231 . Il24l and ||25l (for 
dififerent types of fractional differential equations), which leads to an expression of the solution in 
terms of convolutions of known distributions. In particular we will resort to the class of completely 
asymmetric stable laws (of index less than one). More precisely, let us denote by "p^i-; z), for j - 1,2, 
the density of a stable random variable Xa of index a e (0, 1) and parameters equal to j0 = 1, ;U = 

and cr - {^z\ cos ^) ^ (see If28l for the definitions and the properties of this class of stable laws). As 
well-known, Xa is endowed by the following Laplace transform 

J:{Pa{-,y)\Tl] - , (2.14) 

which will be particularly useful in inverting (12.7b . We need moreover the following result proved in 
II22I : the solution to the following fractional diffusion equation 

g^=c2||, f>0, yeR, ceR 

viy, 0) = 6{y), for < » < 1 • (2.15) 
v,(3;,0) = 0, forl/2<Q'<l 

can be expressed as 

viaiy, f) = — ^ r ^^^^ds = U" {p„(-; m (f), f > 0, y e R, (2.16) 

2cY( \-a)Ja it - s)" 2c 
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where F {■) denotes the Riemann-Liouville fractional integral of order a. By V2a{y, t) we will denote 
the folded solution to (12.15b . i.e. 



2v2a(y, t), y>0 
0, y < 



(2.17) 



Theorem 2.2 The solution to equation ( 12.61 ), under conditions ( 12.51 ), is given, for any k > 0, and 

t>0,by 



Jo 

1 f+cc r 



Pk{y)qv,,v2iy,t)dy 



(2.18) 



I Pv 

Jo 



,(t - s;y)v2v^(y,s)ds + p (t - s;y)v2v2(y, s)ds 



dy. 



k\ JO Uo - JO 

where pk, k> 0, represents the distribution of the standard homogeneous Poisson process N(t), t > 
(with intensity 1) and Py.{-;z) denotes the density of the stable random variable Xy. of index vj G 

■j\y\cos^j 

Proof We observe that ( 12.7b can be rewritten as follows 



rl~v 1 /"I v,-l "2 v,-Upj «1 V, ^ «2 



(2.19) 



Pk(y)' 



dy. 



since, for the distribution of A^, the Laplace transform reads 



■C{pk;j]} = 



1 



-, A:>0. 



(1 +77)*+'' 

The exponential in ( 12.19b coincides with the Laplace transform of the following convolution of the 
stable laws 'p^,^ and p^^: 

8v, ,V2 y) = I Pvii'^-x; y)Py, (x; y)dx. 
Jo 

Therefore, by considering that 



(2.20) 



r(i 



— r 

-V) Jo 



we obtain 



Kit) = 



^r(i - vi) 

«2 

^^r(i-v2) 

»i 

^r(i - vi) 

"2 



wr'j:-' 



Pk(y)e 



/'"- 

Jo 
Jo 

j|(f-w)-"'(^ Pk(y)gy„y2(yv;y)dyyw + (2.21) 

I" -"'"'11 



dy; w I dw + 

X+oo 
Pk{y)e-'-°-^'''''^'^''^'^ydy;w\dw 



Pkiy)gvuv2(w;y)dy\dw. 



Ar(i - V2) 

By inserting (12.20b into ( 12.21b and changing the integration's order, we get 



^r(l-yi) 

'^r(i-v2) 



I Pi'iy^i^J^^ P^S^'^y)^^ J^i^ ~ "^Py^iw - x;y)dw'^dy 



T X ^"^^^(X P^^^^'^^^^' {py,(-^y)}(t-x)dxyy 
+ y £ Pkiy) i^^\{r^y)r- [py,_{-,y)] (f - x)d)^dy. 
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By considering ( 12.161 1 and ( I2.17l i. for c - A/tij, for j = 1,2, formula ( 12.18b immediately follows. □ 



Remark 2.1 The previous result shows that the solution to (12.61 1 can be expressed as the probability 
distribution of a standard Poisson process A^(f), t > 0, composed with a random time argument with 
transition density qv^viiy^ t), that will be denoted as 7"vi,vt (independent from A^): thus we can write 

y,(t) = Pr{A?(7;,,v,(f)) ^k], k>0,t>0. (2.22) 

It is proved in ||2l that the solution to the fractional equation dl.lb is the density of the composition 
of N(t), t > Q with a random time argument Tvit), whose density is given by V2v(y, t). The properties 
of this process have been extensively analyzed in |3 1: it turns out to be a Cox process, with directing 
measure ^qual to A((0^f]) s Ty{f). We will prove below that an analogous result is valid for the 
process A^/,,V2(0 = A^C/"v,,v2(0) introduced here. Moreover we will check that it is also a renewal 
process. 

We derive now a series expression for the transition density q'y, v,(y, f) of the random time- 
argument 7~v,,v,(f)^ f > 0, which is alternative to the integral one given in Theorem 2.2. 

Theorem 2.3 The density qvi,v2{y, of the random time-argument Tvi^viit), f > can be expressed 
as follows 

qvuv2iy,t) (2.23) 



At^^2^rl\ At^^j ^^-v.,i-v2.-v, ^ At^, j+ At^'2 At^' j ^-^-'-"-"^ ^ At^2 



where 



is the Wright function. 

Proof We recall that the solution to the diffusion equation ( 12.151 ) can be expressed as 

V2a(y, t) = {-^-\ , f > 0, y e R, (2.24) 



2ct" ' \ ct", 
(see lfT4ll . for details). Then we get from (12.18b that 

g'y,,v,(y>0 (2.25) 

We now consider the series representation of the stable law of order a e (0, 1) given in lU (formula 
(6.10), p. 583) and already used (with some corrections), in the fractional context, in | 



p„(x;r,l) = -y(-i) 



— : X ^ ^ sm 



^(r + a)(r+ 1) 



(2.26) 



r=0 

In (12.26b the canonical Feller representation for the stable laws (with null position parameter fi) has 
been used, i.e. 

1 /^+°" 

Pjx; y^O^^J e-">' exp }de, a + \; 

hence we must convert the parameters appearing there into those used here, as follows: 



n^Mcos^ n; 

— ; = "tM- 



a - 






2 


7 = 


- — arctan 


n 






f -- 


JIV, 




COS -J- 
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By taking into account the self-similarity property, the stable densities appearing in ( 12.25b become 



A'/ 



l/v, 



(2.27) 



,^ , r(v,r + 1) 

y(-i) — — - 

n\nj\y\) r=Q 



xA^I"' 



s,m(jiVjr) 



= — > (-1) ^-j -TT sm(7rv,r), 

/■=0 ^ ' 



so that ( 12.251 ) reads 



(2.28) 



«!_ y, ,y-i r(v2^+ 1) /n2 



r=0 



sin {nvir) 



1 



'() (f- s)''2''+is''i 



■*'^^-vi,l-vi 



n\\y\ 
As"^ 



ds 



CO 



,_ir(vir+ 1) /«! Y 



^ V " — y sm(7rvir) I ; — 



/■=() 



= TnU-'^ ^'"(^>-^^)L /,r(-v./ + i-vi) Jo (r-.)v.-^'.v.(/.i) -^^^ 



^I^Z^-i) ""^"^■^)L /,r(-v2/+i-v2) Jo (r-.)-.^^'.v.(.i) ^- 



/■=() ^ ' /=() 



1 r(-vir) 



/If'i ^ r! \ .If 

r=() 



r 



«2iyi y 



Uf"-/ ■£;f^!r(l-V2;-V2-Vir)\ /If 2 / 



^Z!r(l-yiZ-vi-V2r)\ if/ 



»2 y j_ /_^^y y 
^fv, Zi H i if. j ^ 



1 



«2lyiy 



where, in the last step, we have used the reflection formula of the Gamma function. □ 

Remark 2.2 Consider the special case n\ - 0, n2 = 1: the distribution of the random order v 
reduces, in this case, to 

n(v) = 6(v - V2), < V2 < 1, (2.29) 

so that equation ( 12.6b becomes the fractional equation (II. lb . with v = V2. The Laplace transform 
(12.71 ) simplifies to 

thus giving the well-known distribution 

pl(t)^A'f^Et:l_,,, (-An, 

for any ^ > (see U for details). The result of Theorem 2.2 can be specialized as follows, for 
ni — 0, n2 — I: 



Jo 



Pk{y)qvuvo{y^f)dy 



-f 

k\ Jo 



/e ^'v2v2{y,t)dy. 



(2.31) 
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In ( 12.3 11 1 we have taken into account that the density of the stable random variable A'y, with /i = 

and 0" = (^ycos ^) degenerates to the Dirac's delta function (i.e. ^y^(w - x;y) - 6(x - w)), 
so that the density ( 12.20b becomes gy^(w;y) - 'Py^{w\y) and (12.211 1 easily yields ( 12.311) . The latter 
coincides with the result proved in |2 1 and already recalled in Remark 2.1. 

As far as Theorem 2.3 is concerned, by putting tii - 0, n2 - 1, the density of the random 
time-argument can be expressed as follows: 

^^.(y, = -^^-y.i-y^ l--^] = ^ivAy^ t), (2.32) 



since in ( 12.231 1 only the term r = of the sum survives. To sum up, the FPP analyzed in (2] is equal in 
distribution to the random time process NiTyit)), whose density can be expressed as a simple Wright 
function; on the other hand, in the distributed order case, the situation is more complicated. The 
density of the process A^(7~y| ,v2(0) is written in terms of infinite sums of Wright functions. Moreover, 
in the single-order case the density V2v2{y,t) coincides with the folded solution to the fractional 
diffusion equation ( 12.15b : in the double-order case the relationship between the density ^vi,v,(3', 
and the fractional diffusion equation of distributed order is more complicated, as we will prove in 
the next section. 

Let us now focus on the probability generating function of the process A/V^vj^ which can be 
expressed in terms of GML functions (11.6b , as the following result shows. 

Theorem 2.4 The probability generating function Gvi,v2(u, t) of the process Nvf,v2 '■^ equal to 

oo 

Gv,,v,(M,f) = Xi«'?I« (2.33) 



7/-+ 1 



^(1 -uy 



+ 



-n- 



m ) ''''2.(>'2-v>)(^+i)+i \ m 

Proof The Laplace transform of Gv,,v2 can be written, by taking into account formula (12.7b . as 
£Ig (u )■ ] - V' " ^ "1^ ^ -yi u A n2T] 



niTf''^^ n2Tf-^^ 



A(\ - u) + niT]^i + n2ri^^ A{1 - u) + niT/"" + n2T]^ 
+ 



(2.34) 



«2 ;|(1 -U)+ |77>'' + TJ-^- -U)+ ^T,". + 77-2 

By applying formula (26) and (24) of [27 1 to the first and second terms of (12.34b respectively and 
recalhng that V2 > vi, the Laplace transform can be inverted as follows: 

«l.v,-v, V^/ "If'-^^'V. .+ 1 / ^(l-«)f' 



y J 1 „ I ^V2,(V2-Vi)(r+l)+l 1 



r=0 



if'-^'Wr^i I A(i-ur^ 



n2 I ''2.(''2-''i)'-+M „2 



which is equal to (12.33b . 
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Remark 2.3 We observe that the infinite sum of GML functions in (12.331 1 coincides with the Fourier 
transform of the solution of the diffusion equation (fractional in time and space), analyzed in IIZTl . 
i.e. 



d 



n, 



^ f(x,t) + '—^f(x,t) ^ c^-^f(x,t), xeR,t>0, 



(2.35) 



in the special case where y = 0. In this case, for - A{\ - u)/n2, equation ( 12.35b reduces to the 
time-fractional equation 



(2.36) 



(with initial condition G{u,Q) = 1). Indeed the probability generating function Gv^^v^_ must solve 
equation (12.36b . as the following steps easily show: 



n2^Gv,,v2(M, + n\—Gy^,y,_{u, t) 



«2^?I(f) + ni^?I(f) 



= ^« 

/t=0 

o 



k=0 k=\ 

-A(\ - u)Gv^^v^{u, t). 



Remark 2.4 By means of the probability generating function, we can check that the distribution 
P^it), sums up to one, for k .... For u - \, formula (12.33b yields 

oo 



k=0 



y i y 



V "2 / r((v2 - vi)r + 1) ^ \ 112 } r((v2 - vi)(r + 1) + 1) 



1, 



since only the term j = in the expression ( 11.6b of the GML function survives. 

Moreover, for u -Q, formula (12.33b gives the probability /5^(0 (akeady obtained in (12.13b ): 



Gv„v,(«,f)|„^o= ^2.38) 



„2 / -y 112 I ^\ "2 / ^"^.("^--'X^+D+l y „2 

/■=0 ^ ' ^ ' /■=() ^ ' ^ 

We note moreover that, in the special case iii - Q,n2 - I, the probability generating function reduces 
to 

Gv,(M,0 = £v,,i(-/l(l 

which coincides with the one obtained for the fractional Poisson process in ||2l, as expected. 

_ We make use of Theorem 2.4 also in the evaluation of the exponential moments of the process 
Nvi,v2 ™d in its resulting characterization as a Cox process. 

Theorem 2.5 The factorial moments of the process 7Vv,,v2, with distribution 'p[{f) and probability 
generating function Gv^^viiu, t) given in ( 12.331 ), are equal to 

E [Ny,,n(t) {Nn,yM - l) ■■■(A?'v,,v,(f) -k+l)] (2.39) 
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Moreover A/Vi.v, is a Cox process with directing measure A ((0, f]) = 7~v^^y^_(t), endowed with density 
Proof We take the k-th derivatives of Gv,,v, with respect to u: 

d" ~ 

—-rGvi,v^{u, t) 



~ h\ "2 ^i!r(v27 + (v2-vi)r+l) (J-k)\ 

«2 j r\^J\r(v2j + (v2-yi)(r+l) + l) (j-k)\ 



which, for u - 1, becomes 

d" ~ 

—-rGvf,v.(u, t) 



u=\ ,•=() 



n2 / r(v2A: + (v2 - vi)r + 1) 
-1 1 + 



"2 / r! r(v2^ + (V2 - vi)(r + 1) + 1) 



-Zl 

Formula (12.40b can be written as ( I2.39K by muhiplying and dividing for k\. 

In order to prove that A/'v,,v2(0> f > is a Cox process with directing measure equal to A ((0, f]) = 
T'vi,v2if), we adopt the characterization of Cox processes by its factorial moments. Indeed it is proved 
in fT2l that for a Cox process they must coincide withjhe ordinary moments of its directing measure. 
Our goal is to show that this equivalence holds for 7Vv,,v, and for the density qy^^^^iy, t) of its time 
argument, i.e. that 

E[T;„v,(f)] = y'qy,,y,{y,t)dy (2.41) 

coincides with ( |2.39t . We start by taking the Laplace transform of ( 12.231 1. which reads 

_ »i 1^ {-n2\y\IX)'' {-ny\y\IX}' T{\ - vyl - V2r - vQ 

-L[qy,.yAy,-),m - ^ Zj rl ljnY(-vil+l-V2r-vi) 7;i-v.'-v2.-v, 

12 ^ (-ni\y\/Ar ^ (-n2\y\/A)' r(l - yj/ - vjr - vz) 



_^«2 y {-ni\y\/A) y 

r=0 1=0 

n^g-mri'Uy\/A - {-n2T]''-\y\/ Af nze"""''"'^''^'* V (-niT]^'\y\/Ay 



nr(-V2l + 1 - Vir - V2) j^l-V2l-v,r-V2 



Z (_-n2T]--\y\/A) «2g y 



so that the Laplace transform of (12.4 U becomes 



(2.42) 

X{E[7;„,,(-)f ;/?} (2.43) 

niA''Tf'-'^k\ n2A''rf^-'^k\ 



(niT]"' +«277*'2)*^' (ni77^i + n2?7''-)*^' ' 
We now take the Laplace transform of (12.391 1. by applying (12.12b : 



/t! : -r + -rk\ 



r,k+l ' I „ \k+\ „k ' / „ \k+\ ' 

"2 (/y^i-"! +f\ "2 (rj"!-"^ + f\ 
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It is simply verified that the last expression coincides with ( 12.431 1. □ 
Remark 2.5 For = 1 we get from (|2.39t the expected value of A/V,,v,: 

EyV„,v,(f) (2.44) 



n- 



n\ jqfl "2 / r(v2./ - Vi7 + 2v2 - vi + 1) 



ni \ «2 / r(v2; - vj + V2 + 1) 



"2-^1 «2 / r((V2-Vi)/ + V2 + 1) 



7=0 

OO 

El 

1=1 

"■2 ^\ «2 / r((V2 - Vi); + V2 + 1) 

Y / mf'--^' y' 1 

ni "2 / r((v2-vi); + V2 + 1) 

/If" 



+ 



"-^y2-vi,y2+l 



"2 " \ "2 

Now consider again the particular case n\ - Q, n2 - l \ formula ( 12.391 1 reduces, in this case, to 



\[NyM{N.,{t)-\)...(Nv,it)-k+\)] 



r(v2/t + 1) 



which coincides with the factorial moments of the FPP obtained in Q. Analogously the expected 
value given in ( 12.44b reduces (for ni = 0, n2 = 1) to 

E7Vv,(f) = — -, (2.45) 

r(v2 + l) 

as expected. We observe that, in the distributed order case analyzed here, both the factorial moments 
and the expected value of A/Vi.vj expressed in terms of Mittag-Leffler functions; for ^ = 1 it 
is a two-parameter Mittag-Leffler function, while, for ^ > 1 , we need a GML function with third 
parameter equal to 1 . This is analogously true, in view of Theorem 2.5, for the ^-th order moments 
of the time argument 7"vi,vi(0- 

We concentrate now on the renewal property of Nvi,v,'- more precisely, we obtain the density 
of the waiting-time of the ^-th event //(f) (or, more exactly, its Laplace transform) and that of the 
interarrival times flit). The latter is expressed again by means of infinite sums of GML functions. 
The same is true for the survival probability *f"'(f). We remark that //(f) can be expressed as the A:-th 
convolution of //(f) and this implies that the process A/Vi.v, is a renewal, since the waiting-time of 
the k-\h event Tk = inf |f > : Nvi^v2(t) = ^} is given by the sum of k independent and identically 
distributed interarrival times Uj, j - l,...k. 

Theorem 2.6 The Laplace transform of the density f^(t) - Pi{Tk e dt} of the k-th event waiting-time 
Tk, is equal to 

■i^{fk'ri]^(-, I, k>l- (2.46) 
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The density of the interarrival time U j is equal to f^,for any j — 1,2, ...andean be written as 

i;-<')^^'--z(-^)'cu,,,.,(-^:i). 

Alternatively 

/-*+oo /-*+co 

flit)- fi(s)g,,.,,(s,t)ds= e-'g,^,,^{s,t)ds, (2.48) 
Jo Jo 

where f\ denotes the interarrival-time density of the Poisson process N (i.e. fi(t) — e^') and gv,,v2 

given in \2.2(M . Then Nvi^v^ is a renewal process with renewal function given by 

»/(f) = £v,-v,,v,+i {-- . (2.49) 



"2 ' \ «2 

The survival probability *!"'(?) = Pr { t/i > f) can be expressed as 



*P''(f)=l- r fiis)ds (2.50) 
Jo 

which solves the relaxation equation of distributed order 

«i^^(o + «2^^(o = -^m (2.51) 

{with initial condition *^'(0) = 1). 

Proof Formula (12.461 1 easily follows from the following relationship 



e-'^'Vx[N,,,,Xt)^k\dt 

/-*+oo 

e-'^'[Pr{n<t}-Pr{n^i<t}]dt 

Jo 



1 

used together with (12.7b . The Laplace transform of the density of the first interarrival time Ui is 
equal to ( 12.46b for = 1 : 

£ {7-; = , (2.52) 

and thus the density of the A:-th event waiting time f1^ is expressed as the fc-fold convolution of f^. 
This proves that A/Vi.v, is a renewal process; its renewal function has been already calculated in 
(12.441) . As a check we show that the well-known relationship between the Laplace transforms of 
m^it) and /j^ holds in this case: 

A 



77 [niT/"' + n2?/''-] 

X{/7;/7} 



A+n\rfi +n2rfi 



The Laplace transform (12.52b can be inverted by applying formula (27) of [271, for a - V2, 
a = ni/n2, P - v\ and b - Ajni, thus giving ( 12.47b . We can rewrite moreover (12.52b as follows: 

£ 1/7; 77] = (2.53) 
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where fi(t), ? > 0, is again the density of the interarrival times for the Poisson process; hence 

X+oo 
/i(fK(T'^-+^'^"^)'^,. (2.54) 

By inverting the Laplace transform (I2.54K taking into account ( 12.211 ) and ( 12.201 ). we get ( 12.481 ). 
We take the derivative of (12.501 ) and we show that -^*I"'(f) = f^{t) given in ( 12.47b : indeed 

at 



^[ m j rl jir(v2j + (V2 - vi)r + I) ^ 



"2 



+ 



fq^l «2 / r! jl^;!r(y2; + (y2-vi)(r+l) + l) 

- fZ.^ „2 j ^''-(>'-''>)'-i «2j «2 j ^--(^-^■x-l)^ „2j- 

The latter expression can be shown to coincide with (12.47b . 

By noting that ( 12.50b is equal to ( 12.331) for m = 0, it is immediately proved that ^'^ solves equation 
( 12.361 ) for u - 0, i.e. equation ( 12.511 ). Alternatively, it is easy to check that the Laplace transform of 
( 12.501 ) is given by 

niT/"'"' + n2'7'''"' 



which coincides with the solution to the Laplace transform of equation ( 12.51b . 



Remark 2.6 In the special case ni = 0, ^2 = 1, from ( 12.47b we retrieve the density of the interarrival 
times of the fractional Poisson process (see [|2|): 

f;(t)^At'^-'Ey,,y,(-At'''-). (2.55) 

Likewise the survival probability (12.50b reduces to 

»I"'(f) = £v„i(-^f''^). 



It is interesting to analyze the asymptotic behavior of the waiting time densities and of the re- 
newal function and to compare these expressions with the corresponding formulas obtained for the 
fractional Poisson process. To this purpose we need to prove the following integral representation 
for the GML function: 

E\A,-cf) = r^'-l^e-" ' ^ ' ' (r + ce ) ^^^^^ 

27ri Jo [r^" + 2r''ccos(v7T) + c^] 

We start by checking that, for k = I, formula (12.56b coincides with the form given for Ev^i—t^) in 
11, i.e. 



E.,p{-ct') = r'-^e-" \ ;\ \ ' dr. (2.57) 

" ' r" + Ir^c cos(v7!-) + 



14 



In order to prove formula (12. 56) we multiply and divide the m-th term in the series expression of 
E'^pi-cf) for sin((j0 + vm)n)lji and apply again the reflection formula of the Gamma function, as 



follows 



(2.58) 



1 -^{m + k- \)\{-cfT sin((j6 + vm)n) 



m\Y{vm + B) 

m=\J 

Tik-l) 



r(l - vm -/3)r(vm+/3) 



n{k-\)\^^ m\Y(vm+l3) ''X J,, 

n{k-\)\]^ Jo {^^ m\T{vm+li) 



(m + k- ly.i-cy")'" e'"™^'"!^ - e-'"""-'"!^ 



dy 



2m(k- 1) 
= [by ([231] 

= — e-" 
Jo 



Xoo /-»oo 



dr. 



This coincides with ( 12.561 ). The asymptotic behavior of E'^^^i-cf) can be obtained from ( 12.561 1 and 
reads, for f — > oo: 



-vk /-oo 



f 



f 



+ 2c^ cos TTV + C^) 



2\* 



I^i-^sin(.(/^-v.))..(r-) 



1 



(2.59) 



c*f''*r(y8 - vk) 



f ^ oo. 



For k = I, formula (12.591 1 reduces to the one holding for the Mittag-Leffler function, which can be 
deduced by (12.57b . i.e. 



1 



ctTiP - v) 



+ 0(0, t 



(2.60) 



For f — > 0, we get instead that 



1 r°° 

lim E^gi-cf) = lim — e"'; 



(z^v + 2czT cos;ry + c^fi^y 



dz 



^ Jo 



2v* „sm(;i^ r(l-^) . 1 
e ~z ^ t/z = sm(;rje) = 



r(j8) 



(2.61) 



The asymptotic behavior for small f can be deduced directly by the series expression of E'^p{-cf): 
indeed we get, for < f « 1, 



1 _ cfk 



which reduces, for A; = 1, to the well-known expression (see ll20l . formula (3.13)). 
The interarrival-time density (12.47) can be rewritten, by applying ( 12.581 ). as 
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"2 



(2.62) 



CO I 

r=0 ^ 



"2 



27ri 



. I 

Jo 



-z'2V2(r+l)-V2-(v2-vi)r 



g/;rV2+/;r(v2-vi)r ^-CTV2-/;r(v2-vi)r 



,/;ry2 



e ~'dz 



2m 



w — 



Jo nie''r(i'2-vi)2Vi + „2Z''2 + Ae"^2 



Jo nie-"< 



e ''dz 



Jo 



„je'>r(V2-Vl)^^V,fl-V| +„2W^2fl-V2 +^e'>V2f 



n2r(v2) ' 



which shows that, for f — > 0, the asymptotic behavior of depends only on the larger fractional 
index V2- The same conclusion can be drawn by looking at the series expansion of /j''(f) given in 
(IT471) . 



For t 



, from the sixth line of (12.621 1. we have instead that 

m 



(2.63) 



2ni 



e-"'dw 



(„jeiMv2-v,)^v,f-v, +„2W^2r''2 + Ae''^^2)(«ie-'''(''2-''i)w^ir''' + n2W''^t-''^ + Ae-'""'-) 
nivir'""' sin(;rvi)r(vi) _ niVir'""' 



An 



f ^ 00, 



Ar(l-yi)' 

which depends only on the smaller fractional index vi Both asymptotic expressions (I2.62l i and ( 12.631 1 
are exactly the same as for a fractional Poisson process of single order equal to V2 and vi respectively 
(see 0, formulae (2.38) and (2.36)). 

Analogously we can analyze the asymptotics of the survival probability, which turns out to be, 

f ^ CO, 



for 



2m Jo nie-''^(^2-^i)w^ir''i +n2w''-r^^ 



(2.64) 



+ 



„^j-v,g;;r+/;r(v2-v,) 

27r/ 



e-^w^'-^dw 



27r/ Jo «ie-''^(''2-*'i)w''ir''i H-nivV^f 

niT''' sin(7r(l - vi))r(vi) _ niT" 

For f ^ 0, by writing down the first 



12 w''2r''2 + Ae-'""'- 
f — > 00. 



Ar(l-yi)' 

terms of the series expansion in (12.501 ) (at least for j = 0,1,2 
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and r = 0, 1 , 2) and doing some manipulations, we finally get 

TO - +0(1"'), 0<f«l. (2.65) 

«2 r(V2 + 1) 

As far as the renewal function is concerned, its asymptotic behavior can be represented as follows: 

m\t) - — -, t ^ oo (2.66) 

niHvi + 1) 

and 

m\t) - — -, t ^ 0. (2.67) 

n2r(v2 + 1) 

From ( 12.661 1 it is evident that the mean waiting time, which coincides with lim,^ootlm'{t), is infinite, 
since vi < 1. 

Remark 2.7 We remark that, in Theorem 2.6, the survival probability *P'' expressed in (12.50b is 
proved to be a solution of the relaxation equation of distributed order (12.51b under the double-order 
hypothesis (12.5b . This result can be compared to the analysis in lfT6l . where only the Laplace trans- 
form of the solution is presented, together with its asymptotic behavior. Formulae (12.64b and (12.65b 
coincide with the result (4. 13) obtained therein, but here we provide an explicit formula of the solu- 
tion, in terms of infinite sums of GML functions. 

2.2 Interpolation between fractional and integer-order equation 

We analyze now the following equation: 

«i^+«2% = -^(/?*-PA-i), ^>0, ve(0,l) (2.68) 
df dt 

which is obtained from ( 12.6b . as a special case for V2 = 1, under the usual initial conditions 

pm-[l ll\ , (2.69) 

and p-\{t) - 0. Equation ( 12.68b represents an interpolation between the standard and the fractional 
equation governing the Poisson process. Hence the solution, which will be denoted in this case by 
p^, must coincide, for n\ = 0, «2 = 1 with the distribution of the homogeneous Poisson process, i.e. 
Pk,k > 0. On the other hand, for n\ - 1, n2 = we must retrieve the distribution of the fractional 
Poisson process, i.e. pJJ, k>Q), given in (11.5b . 

The Laplace transform of the solution to equation (12.68b can be obtained directly by putting 
vi = V and V2 = 1 in the result of Theorem 2. 1, so that we get 

■L\PL,ri\ — -, k>0. (2.70) 

for any A: > 0. In order to invert this expression we adapt the result of Theorem 2.2 as follows. 

Theorem 2.7 The solution to equation \2.68h under conditions \2.6% , are given, for any k > 0, and 
t>0,by 

r*+oo 

Pl(t)^ Pk(y)qAy,t)dy (2.71) 
Jo 

" ^ X " + n2P,{t-y))] dy. 

Here "Pyi-^y) denotes the stable law of the random variable Xy of index v € (0, 1) and parameters 
equal to /3 — I, fi — nxly]/ A and cr — (^lyl cos 
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Proof We observe that ( I2.70l i can be rewritten as follows 
so that we get, by an argument analugous to Theorem 2.2, 

+^ Pk(y)gv(t;y)dy 
^ Jo 

The density in (12.721) can be expressed as follows 

pjw - x;y)6{x —)dx 

Jo ^ 

= Pv(w;yy, 



hence 



Jo 

which coincides with ( 12.711 ). 



—V I n\ i _v= "2 = 

PaXO = Pkiy) -T^, r (t-w) p^(w;y)dw+ —p^{w;y) 

Jo L^^ -V) Jo A. 

dy. 



dy 



^r(Py(w-y)) + ^Py(w;y) 

A A 



Remark 2.8 For 112 -Q and ni - 1, we get that Py{t;y) = Py(t,y). Therefore formula (12.71) reduces 
to 

?I(0 = ^J^ /e-yripy(-,y))it)dy 

y'e-yv2y(-,ym)dy (2.73) 



^ Jo 



as for the single-order fractional equation (see (12. 3U and ( I2.16I )-( ITT7] )). On the other hand, for 
111 - and 112 - 1, it is Py(t',y) = S(y) and /7^(f) = pa:(0. since, in this case, equation (12.681 1 reduces 
to the equation governing the Poisson distribution. 

Remark 2.9 A particular feature in this section is that for the process governed by ( 12.681 ) the prob- 
abiUty generating function Gy, as well as the probability of zero events /JJJ and the interarrival time 
density f^(t), can be expressed as infinite sums of the Kummer confluent hypergeometric function 
iFi ia,j3; x) . The latter is defined as 

i^i(«;r;z) = 27T^' ^^^^c, rec\Zo, (2.74) 

where (y)^ = y{y + l)...{y + r - I) (for r - 1,2, and y 5^ 0) and (y)o = 1, or, in integral form, as 
iFi(g;y;z)^ ^^/^ C f-'il - t)y-"-' e'' dt, < 5R(a) < 5R(r), (2.75) 

r(Q')r(r - a) Jo 

(see ||9l, p. 1085). Indeed it is well-known the following relationship between the GML function with 
first parameter equal to one and i^i (a, y; x): 

EiJD^-^^iFi ia;y;z) 
r(y) 
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(see E], p.62). 

By specializing result (12.331 1. the probability generating function Gv(u, t) is equal to 

■^r(r(l-y)+l)' H «2 / 

Zjr((r+l)(l-v) + l) ' H «2 
Analogously, from ( I2.47l i the interarrival time density reads, in this case, 

= - Z rr'^^L h 1; ^ - + 1; . (2.76) 
"2 fqf r(r- vr+ 1) \ n2/ 

Remark 2.10 The expected value of the renewal process Nv(t), f > with distribution (12.711 1 is given 

by 

UNyit) = -£i-.,2 (-^^i^) , (2.77) 

"2 \ "2 / 

SO that we get this asymptotic behavior 



EM(f) - — — -, t ^ oo. 

n\Y[l + v) 

This expression coincides with (12.451 1. for n\ - \: the mean value is not influenced by the presence 
of the first derivative. On the contrary, for f — > 0, we obtain from (12.77b that 

E7Vv(0 ^ — , f ^ 0, 

"2 

i.e. the usual expected value of the Poisson process. Therefore the first derivative dominates equation 
(12.681 1 asymptotically, as f — > 0. 



3 Diffusion equations of distributed order 

We study equation (11.13b in the double-order hypothesis (12.5b . i.e. 

«i— - + n2^— = 7-7, X € R,f > 0, v(x,0) = 5(x), ni,n2 > 0, (3.1) 

for < vi < V2 < 1 . Equation (13.1b can be viewed also as the particular case (for y = 2) of equation 
( 12.35b analyzed in 1,27 J . In that paper only the Fourier transform of the solution is given in explicit 
form, in terms of infinite sums of GML functions. Our aim is to give an explicit form of the solution, 
by using an approach similar to the previous section and providing an expression of the density of 
the random time in the subordinating relationship (I1.14t . This turns out to coincide with the density 
of the random time Ty^^v^{t), i.e. with qv^.vi given in (12.18b or (12.23b . 



Theorem 3.1 The solution to equation ( 13.71 ). is given by 

f(x,y)qvi,v2(y^t)dy (3.2) 



dy. 



I . I Py^t - s;y)v2yfy,s)ds + I 

Jo -JATTy [Jo ' Jo 



Pv. (t-s; y)v2v, {y, s)ds + Py,(t - s; y)v2v, (y, s)ds 
Jo 

where f is the transition density of a Brownian motion B(t), f > 0, V2v, is given in A2.16\l - ^.17]l and 
Py.(-'-,y) denotes the stable law of the random variable Xy. of index vj e (0, 1) with parameters /3 — I, 
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jU = and cr — (nj\y\ cos ^ ^ ,/or j — 1,2. Alternatively the density qv^.Vl in ( 15.21 ) can /je written 
as in Ii2.23\l . 

Proof We take the Fourier transform of (13. U . so that we get 



y(6>,0) = 1 



(3.3) 



where ^ 

Vy,,yMt)^r{Vy„y,;0} = r 

Taking now the Laplace transform of (13.3b we get 



X{yv,,v,;0,77} 



(3.4) 



We can invert the Laplace transform, by noting that it coincides with (12.1 11 1 for A - 6^, as follows: 



K„v,(0,f) 



(3.5) 



- z- 



nit" 



"2 



^l'2,(V2-Vi)r+l 



"2 



^V2,(V2-Vl)(r+l)+l 



g|2jV2N 
"2 ; 



thus obtaining a first form for the solution to (13. 3t . Since inverting the Fourier transform (13. 5t seems 
not possible in closed form, we rewrite (13.4) as follows: 



1 r+°° 



We note that the term e ("ii'^+nrn'^^ ^.^^ ^^^^ again as the convolution of two stable laws Py, of 
index vj e (0, 1) and parameters equal to yS = I, /j. - and cr - Qnj|w| cos ^) ^ ^ for j - 1,2 (see 



(12.14b ). Therefore we get, alternatively to (13.5b 

Vy„y,(e,t) 

r (^-z)-* r 

r(i-vi)Jo Jo 

Jo Jo 



Jo 



(3.6) 



Py{z-x;w)p{x;w)dx 



r(l-V2) 

-•+00 



(z-.x;w)/?^,,(x;w)t/.x 



dw + 
dw 



X+oo r w 

e-"'"' r' {Fv, (■; vf)} (^) :Pv2(f - ^; v^)^?^ 

X+(X5 r 



dw + 



(t — x; w)dx 



dw 



V2v,(x, w) p (t- x; w)dx 



dw + 



Jo [Jo 



V2v2(x, w) py (t - x;w)dx 



dw. 



where again V2v, is the solution to equation (12.15b with = ^/nj, j = 1,2. Finally, we recognize 
in ( 13.6b the Fourier transform of the Gaussian density, with variance 2\w\, so that we can write the 
subordination relationship (13.2b . □ 

The previous theorem shows that the solution to the double-order equation (13.1b can be seen as 
the density of the random-time process 



(3.7) 



20 



where B is a Brownian motion (with infinitesimal variance equal to 2) and 7~v,,vt is the random time, 
independent from B, with density ^y, y, given in (13.21 ) or, alternatively, in (12.231 1. for A - 1. By using 
the results obtained in Theorem 2.5, we can evaluate the moments of Sv,,v2, as follows: 

Theorem 3.2 The k-th order moments of the process ^v^yi '^^^ given by 



Vi,V2 

CO for ;t = 2/1 + 1 



""+' ^ V2-V1 ,V2"+V2-Vi +1 V Hi / nJ, V2-Vl,V2/J + l \ "2 / 



for k = 2h 



Proof By the definition < l3.7l t we can write 



ES* „^(f) = E[B(T;,,y,(f))] = J f(x,y)qy,,y,(y,t)dydx. 

The odd order moments of ^^, ^2 are obviously null, while the moments of order 2h, for h eN, can 
be evaluated as follows: 



o2h 



(f) 



r 

?V,,V2(3'>f) 



= — ^vi.v^C}',?) (43;w)'' — t/wt/y 
Jo Jo yw 

22h I Y\ r+°° 

= -^r \h + - j /?v,,v,(3', f)'^}' 

r*+oo 

y''qv,.v,(y,t)dy 

Jo 



= 2 



(2/1-1)! 

(h-iy. 



-{2h)\El^l_ 



V2-V1 ,V2/i+V2— Vi + 1 



+ ^(2/i)!£,t', 



■2/1+1 



"2 



where, in the last step, we have applied formula (12.391 1 and the relationship 

E [T.,.yM^ = E [Ny„yM {Ny„y,(t) - l) ...(Ny„y,(t) ' k + I)] , € ! 

proved in Theorem 2.5. 



(3.8) 
□ 



Remark 3.1 We can check the previous result by noting that, for h = 1, we get the second-order 
moment obtained in [jSj (see formula (16), for t = D = 1): 



2«if2^^-''' 



(3.9) 



(; + l)!(-^y 2^-". ^ (;>l)!(-^y 



"2 ^ ;!r((v2-vi)7 + 2v2-vi + 1) n2 j^^ jlT ((v2 - vi) j + V2 + I) 



^ ;!T 



2?" 



2?" 



"2 -f^ r((v2 - Vi)/ + V2 + 1) n2 ^ r((v2 - Vi)y + V2 + 1) 



2f^ 



"2 



y2-vi,V2+l 



"2 
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Our attention is now addressed to the equation solved by the density ^v,,vt of the time argument 
7~y|,v, (which is shared by the processes TV^i vj andSv,_v2 )■ In analogy with the single-order fractional 
case this equation must be of "second-order" (involving the two fractional indexes vi , V2), but is not 
evidently given by ( 13.11 ). We prove in the next theorem that a further time-fractional derivative must 
be included in the diffusion equation (13. Il l in order to obtain ^v,,v, as solution. 

Theorem 3.3 The density qy^^y^(x, i) coincides with the folded solution 

i 2v(x,t), x>0 



of the following equation 



"'5^ ^"'5^^ 5?' ^^^'t>Q,nun2>Q, (3.11) 



for < Vi < V2 < 1, with initial conditions 

v(x, 0) = 6(x), for < vi < V2 < 1 
|v(jc, f)|,^o = for i < vi < V2 < 

Proof Let us take the Fourier transform of (12.42) . which reads 



(3.12) 



e"'J:{qy,.y,;7]}dy 

CO 

- («l7?^'-'+«2/7''^-') 



(3.13) 



1 



dy 



1 



nii]^' + n2ri^^ - iO niT]^' + nirj^- + id 
2(«i77''' +«2'7''^)2 



by taking, for simplicity, A - \. We take now the Laplace-Fourier transform of equation (13.111 1. by 
considering formula (12.8b together with the initial conditions ( 13.121 ): 



+2n,n2rf^^'^J:[Qy^,y^■e,r^}-2nln2lf^^''--' = -fl'X {ev,.v,; 0, /?) , 
whose solution coincides with (13.131) . by taking into account ( 13.10b . □ 

Remark 3.2 Equation ( 13.1 lb can be itself interpreted as a distributed order fractional equation, by 
assuming a different density for the random fractional index v, which, in this case, is defined on the 
interval (0, 2]: indeed we can formulate n(v), as follows: 



n(v) — n^6(v - 2vi) + n^6(v - 2v2) + 2nin26(v - (vi + V2)), 



< Vi < V2 < 1 , 



(3.14) 



for ni,n2 > and such that ni + n2 - 1. The last condition is enough to fulfill the normalization 
requirement Wj + + 2ni«2 - 1. 

By considering ( 12. 39b together with (13.8b we can obtain the second-order moment of the diffusion 
process Tv,_y, governed by equation ( 13.1 lb . i.e. 
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E [r,,,(0] = ^;^£tv„3v.-v,.l (-^) + ^£tv„2v... (-^) (3.15) 
«lf3v2-v, - (j + 2)(j+l) 



z 



Wo r((v2 - vi); + 3v2 - vi + 1) \ n2 

^ ^ 0- + 2)0-+l) / wif^-^'-V' 



-— y 

«2 ^ r((v2 - vi); + 2y2 + 1) \ n2 
[' = 7 + 1] 

«2 it r((v2 - vi)/ + 2v2 + 1) i 112 )^ nlj^ 



(/ + 2)(Z+ 1) 



r((v2 - Vi)Z + 2V2 + 1) \ «2 



.f2 I 



2 vt-V|,2vt + 1 I 

\ "2 

We compare ( 13.151 1 with the second-order moment of the diffusion process Sv,_y, governed by equa- 
tion (13. lb (which is given in ( 13. 9t ). noting that, apart from the similar structure, the second-order 
moment of Tvi,v2 is expressed in terms of a GML function, instead of a standard one. For f — > oo the 
different asymptotic behavior is described below, by applying (12.591 1 and (12.601 1: 

^ J "^(1 -I- 2vi) 



?v,.v,(0 



If 



mm +vi) 



For f — > 0, since the behavior of the GML and the standard Mittag-LefHer function coincides, we 
must apply, in both cases, formula (12.611 1. Hence we get: 

L J n^r(l + 2v2) 



n2T{\ + V2) 



We can conclude from ( 13.16b and (13.17b that, in the case where v\,V2 < 1 /2, the effect of diffusion 
with retardation, which is characteristic of Sv,,v, (see also 161), is emphasized for the process 7"v, _v,, 
since the difference between 2v2 and 2vi is greater than for V2 and v\ . 

A different conclusion should be drawn in the case where either V2 or both v\,V2 are greater 

than 1/2. For v\ < 1 12 and V2 > 1/2, we can observe that the asymptotic behavior of E [Tv,,v2(f)] 
drastically changes, at least for f — > 0: indeed in this case it goes to zero faster than t. For vi , V2 > 1/2, 
in addition to this effect, we note that, also for f ^ 00, the rate convergence is greater than in the 
standard diffusion case (besides that of diffusion with retardation). Therefore the process governed 
by (13.1 lb . for vi, V2 > 1 /2, can be interpreted as a diffusion with acceleration. 
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